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The NA38 and NA50 experiments at the CERN SPS have measured charmonium 
production in different colliding systems with the aim of observing a phase tran- 
sition from ordinary hadronic matter towards a state in which quarks and gluons 
are deconfined (Quark Gluon Plasma, QGP). In fact it was predicted that the 
J/ip yield should be suppressed in deconfined matter. The analysis of the data 
collected by the NA50 experiment with Pb-Pb collisions at 158 GeV/c per nu- 
cleon shows that the J/ip is anomalously suppressed in central collisions and the 
observed pattern can be considered as a strong indication for QGP production. 
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1 Introduction 



In ordinary nuclear matter, quarks and gluons are confined inside nucleons. 
Non-perturbative calculations of Quantum Chromo Dynamics predict that, 
when temperature exceeds a critical value T c ~ 150-180 MeV, nuclear matter 
should undergo a phase transition into a state of matter in which quarks and 
gluons arc no more confined into hadrons and behave as free particles. Such 
a state of matter is named Quark-Gluon Plasma (QGP). 

Heavy-ion collisions are a powerful experimental tool to investigate nu- 
clear matter under extreme conditions: the formation of the QGP is expected 
to occur in these collisions if the critical temperature and energy density re- 
quired for the phase transition are reached. Several probes have been proposed 
as signatures of the formation of a deconfined state of matter. In particular 
Matsui and SatzEl predicted that the J/ip yield would be suppressed in a de- 
confined medium due to the Debye screening of the attractive colour force 
which binds the c and c quarks together. J/tp suppression is a particularly 
interesting signature of QGP formation because it probes the state of matter 
in the earliest stages of the collision, since cc pairs can only be produced at 
that time. Moreover the J/ip is a tightly bound state that can not be easily 
broken by interactions with the hadronic medium and therefore it carries its 
original message through the different stages of the reacting medium. 

At the CERN SPS, the NA38 and NA50 experiments have studied J/tp 
production using the [i^ [i~ decay channel with incident proton, oxygen, sul- 
phur and lead ions on several targets. 

2 Experimental setup and data taking conditions 

The J/ip is detected via its [V fi~ decay by a dimuon spectrometer which 
consists of an air-gap toroidal magnet equipped with two sets of multiwire 
proportional chambers (4 MWPC upstream the magnet and 4 downstream) 
for muon tracking purposes and with 6 scintillator hodoscopes to provide the 
dimuon trigger. The covered rapidity window is 2.8 < yi a b < 4.0. 

The centrality of the collision is estimated event by event using three dif- 
ferent detectors. The electromagnetic calorimeter measures the neutral trans- 
verse energy (Et) produced in the interaction in the pseudorapidity window 
1-1 < rjiab < 2.3. The zero degree calorimeter measures the energy Ezdc, 
essentially carried by the projectile spectator (non-interacting) nucleons. It 
covers the pseudorapidity interval r\i a \> > 6.3. Finally the silicon microstrip 
detector measures the multiplicity and the angular distribution of charged 
particles in the acceptance window 1.5 < r\i a b < 3.5. 
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A quartz beam hodoscope (BH), placed about 33 m upstream the target 
is also used to count the incident lead ions and to reject beam pile-up. The 
NA50 target is a segmented active Pb target: there can be up to 7 subtargets 
separated by few centimeters of air. Two quartz blades located off the beam 
axis on the left and right side of each subtarget identify where the interaction 
has occurred. For a more detailed description of the detectors see Ref. 4. 

3 Charmonium production in p A and light ion interactions 

NA38 and NA51 experiments collected an extensive set of measurements of 
J/ip production using p, O and S beams on several targets. Since charmonium 
production is a hard process and therefore it scales as AxB (where A and B 
are the projectile and the target mass number respectively), it is possible 
to define a cross-section per nucleon-nucleon collision as B w crj/^,/(vl x B), 
where B w is the branching ratio of J '/ip into two muons, in order to compare 
charmonium production in different colliding systems. 

The standard analysis method has been described in detail in Ref. 4. The 
number of J/ip events is determined by fitting the invariant mass spectrum of 
opposite sign muon pairs above 2.9 GeV/c 2 . The fit is performed including 
5 contributions: J/ip, ip' , Drell-Yan, open charm (i.e. semileptonic decays of 
D and D mesons) and combinatorial background from tt and K decays. An 
example of invariant mass spectrum is presented in fig. [I] 

In figure | the J/ip cross-section per nucleon-nucleon collision as a func- 
tion of AxB is represented: it can be seen that all the data from p-p up 
to central S-U collisions show a continuous and monotonia.. J/ip suppression 
pattern from the lighter to the heavier interacting nucleiod. The systemat- 
ics of J/ip production from p-p to S-U canJpe parametrized by the simple 
law cjab = (TQ{AB) a with a = 0.918 ± 0.0150. Within the framework of the 
Glauber model, using a simple first order exponential fit, these data lead to 
a J/ip absorption cross-sectiorLof 5.9 ± 0.6 mb (equivalent to 6.4 ± 0.8 mb 
for a full Glauber calculation)!! This result can be understood in terms of 
ordinary nuclear absorption of a preresonant cc state meant to become later 
on, if not destroyed, the fully formed J/ip. This normal suppression sets the 
baseline with respect to which we can compare the pattern of J/ip production 
in Pb-Pb interactions. 

4 Anomalous J/ip suppression in Pb-Pb collisions 

In fig. [2] it can be seen that the J/ip cross-section measured in Pb-Pb collisions 
lies ~ 5 standard deviations below the value expected from the fit to the data 
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Figure 1. Opposite sign muon pair invariant mass spectrum for Pb-Pb collisions. 



from p-p to S-U and this result indicates that there is a new suppression 
mechanism at work. 

In order to perform a study of charmonium production as a function of the 
collision centrality, it is necessary to find a replacement for AxB in the defi- 
nition of the cross section per nucleon-nucleon collision. Since the Drell-Yan 
(DY) process provides muons in the same mass range as the J/ip and its cross 
section is expected to be proportional to the number of elementary nucleon- 
nucleon collisions without any sizeable nuclear effect, the ratio B^^aj/^JaDY 
can be used to study the centrality dependence of the J/tp suppression!^ with 
the advantage that in the ratio oj/^/o'dy systematic errors related to detector 
inefficiencies and flux uncertainties cancel. The transverse energy spectrum is 
then divided into bins: in each of these bins a fit to the invariant mass spec- 
trum of the muon pairs is performed and the J/tp over DY ratio is calculatedQ. 

The disadvantage of this analysis method is that there are statistical 
fluctuations in the o~j/^/o-dy ratio, essentially due to the small statistics of 
the DY sample. In order to overcome this problem, a new independent analysis 
method has been developed!!! in which the sample of DY events is replaced by 
the huge sample of minimum bias (MB) events, collected with a beam trigger 
which fires every time that a non zero energy deposit is detected by the Zero 
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Figure 2. J/ifi cross section as a function of AxB from p-p to Pb-Pb interactions. 



Degree Calorimeter (even for the most central collisions some energy, from 
produced particles, falls into the ZDC acceptance). 

In this anaysis method the number of J / ip events is obtained without any 
fitting procedure, simply by counting the number of opposite sign muon pairs 
in the mass range from 2.9 to 3.3 GeV/c 2 , after combinatorial background 
subtraction. In order to compare the results of this analysis with the ones 
obtained with the standard method, the minimum bias reference has to be 
converted into the Drell-Yan reference. This is done multiplying the measured 
Et spectrum of the MB events by the ratio of the theoretical shapes of the 
Et spectra of DY and MB, evaluated in the frame of the Glauber model. 

Because of the bigger size of the minimum bias sample, the new analysis is 
no more affected by statistical fluctuations and a finer Et binning is possible. 
It is also free from most inefficiencies as it is still computed from a ratio of 
experimental numbers. 

The J /ip over DY ratios obtained with the two independent analysis show 
a good agreement, as it can be seen in fig. ^ where the complete J/ip sup- 
pression pattern obtained with the standard (closed circles) and the minimum 
bias (open points) analysis for the 1996 and 1998 data samples is plotted as 
a function of the transverse energy E^. 
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Figure 3. B^^aj a/, /crnY ratio as a function of Et, obtained with the standard and mini- 
mum bias analysis of the 1996 and 1998 data samples. 



The solid line visible in the figure represents the "ordinary" nuclear ab- 
sorption of the cc pair, which accounts for the J/ip yield measured in lighter 
colliding systems (from p-p to S-U). It can be observed that for the most 
peripheral Pb-Pb collisions the J/ip suppression can be accounted for by or- 
dinary nuclear absorption while a clear departure from this trend can be seen 
at Et ~ 40 GeV, which approximately corresponds to an impact parameter of 
8 fen. It is also evident that there is no saturation of J/ip suppression for the 
most central Pb-Pb collisions: a second drop in the J/ip suppression pattern 
can be observed at Et ~ 90 GeV. 

5 Conclusions 

The combined results of the NA38 and NA50 experiments clearly indicate a 
step- wise J/ip suppression pattern with no saturation for the most central 
collisions. 

The experimentally observed suppression is much steeper (when the ex- 
perimental resolution is taken into account) than the results of the pre- 
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dictions of all presently available models of J/ip suppression based on the 
absorption of the meson by interactions with the surrounding hadronic mat- 
ter. 

On the contrary the observed suppression pattern can be naturally un- 
derstood in a dcconfincmcnt scenario. Since w 30-40% of the J ftp's come 
from radiative decay of \c mesons, the first anomalous step in charmonium 
suppression can be understood as due to the melting in dcconfincd matter of 
the Xc which is less tightly bound than the J ftp. In this scenario, the second 
drop, observed for more central collisions, would be due to the dissolving of 
the more tightly bound J/ip meson, which requires energy densities grater 
than the ones necessary to melt the \ c . 

Hence it can be concluded that the J ftp suppression pattern observed in 
the NA50 data provides significant evidence for dcconfincmcnt of quarks and 
gluons in Pb-Pb collisions. 
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